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£^ , Abstract. Massive stars and supernovae (SNe) have a huge impact on their 

■ environment. Despite their importance, a comprehensive knowledge of which 
massive stars produce which SNe is hitherto lacking. We use a Monte Carlo 

— . method to predict the mass-loss rates of massive stars in the Hertzsprung- Russell 

J> ■ Diagram (HRD) covering all phases from the OB main sequence, the unstable 

Q\ | Luminous Blue Variable (LBV) stage, to the final Wolf-Rayet (WR) phase. Al- 

though WR produce their own metals, a strong dependence of the mass-loss rate 

■ on the initial iron abundance is found at sub-solar metallicities (1/10 - 1/100 
solar). This may present a viable mechanism to prevent the loss of angular 
momentum by stellar winds, which could inhibit GRBs occurring at solar metal- 

■ licities - providing a significant boost to the collapsar model. Furthermore, we 
discuss recently reported quasi-sinusoidal modulations in the radio lightcurves 

r"| . of SNe 2001ig and 2003bg. We show that both the sinusoidal behaviour and 

i' the recurrence timescale of these modulations are consistent with the predicted 

' | mass-loss behaviour of LBVs. We discuss potential ramifications for the "Conti" 

scenario for massive star evolution. 

£3 ■ 

c/3 . 
Ctf ■ 

> : 

■ 1. Introduction 

><: 

Massive stars have a huge influence on their environments via stellar winds 
and their final explosions. However, we currently do not know with any degree 
of certainty, which massive stars produce which type of supernova (SN). The 
evolution of a massive star (M > 40 Mq) is generally believed to be driven by 
mass loss, as described in the "Conti" scenario: O — > Luminous Blue Variable 
(LBV) — > Wolf-Rayet (WR) star (Maeder, this meeting), where the WR stars 
include both nitrogen-rich (WN) and carbon-rich (WC) stars. 

Mass loss determines the stellar mass before collapse and is thus relevant 
for the type of compact remnant that is left behind (i.e. neutron star or black 
hole). This process is expected to d epend on the metal content (Z) of the host 
galaxy (e.g. Eldridge Sz Vink 20061 ). Furthermore, as WR stars are the likely 



progenitors of long-duration gamma-ray bursts (GRBs), the strength of WR 
winds as a function of Z is especially relevant for setting the threshold Z for 
forming GRBs. 

Given the crucial role that mass-loss plays for massive star evolu tion, we 



have computed mass-loss rates using a Monte Carlo method, described in lAbbott Sz Lucy 
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(119851 ). Ide Koter et al.1 (119971 ) and I Vink et~afl ([1999ft . We discuss the predictions 
in Sects. 2-4 in order of decreasing temperature: WR stars — > OB supergiants 
— ► LBVs. In Sects. 5 and 6, we link our mass-loss predictions with certain types 
of radio SNe. Conclusions and outlook are in Sect. 7. 



2. Wolf-Rayet mass-loss rates as a function of metal content 

In recent years, it has become clear that gamma-ray bursts (GRBs) are associ- 
ated with the final e xplosion of a massive star, p roviding enormous impetus to 
the collapsar model ( MacFadven fc Wooslev 1999). The model works best if the 



progenitor fulfils two criteria: (i) the absence of a thick hydrogen envelope (so 
that the jet can emerge), and (ii) rapid rotation of the core (so that a disk can 
form). This may point towards a rapidly rotating WR star. 
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Figure 1. Mass loss versus initial Z for late- type WN stars (solid line) and 
WC stars (dashed line). Note that se lf-enrichment is account ed for, but does 
not enter in our expression of Z. See lVink fc de Koteri J2005) for details. 



In the so-called "Conti" scenario (IContilll976h . WR stars are the result of 
mass-loss during earlier evolutionary phases, while in a complementary scenario, 
the removal of the hydrogen envelope may be due to a companion. Recently, 
an alternative scenari o for producing a GRB progenitor has gained popularity 
(|Yoon fc Langeril2005l : IWooslev fc Hegerl 120061 k when a star rotates rapidly, it 
may mix "quasi homogeneously", and evolve along a track that more or less 
coincides with the zero age main sequence. A problem for producing a GRB 
within this scenario however, is that the WR stars in the Galaxy possess strong 
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winds, which may remove the angular momentum ( Langer 19981 ). making it 
challenging, if not impossible, to produce a GRB at Galactic Z. 

This however might not be an issue if WR winds are weaker at low Z, so 
the question is: "are the winds of WR stars Z-dependent?" and if so, "how 
strong is this depe ndence?" The dense w i nds of WR stars are likely driven by 
radiation pressure ( Nugis <fe Lamersl 2002 ; Grafener Hamann 20051 ). just like 
their less extreme O star counterparts, which have been known to be driven by 
radiation pressure since the early 1970s. This in itself need not necessarily imply 
that WR winds depend on metal content. WR stars produce copious amounts 
of metals such as carbon (in WC stars). If, on the one hand, these self-enriched 
elements dominate the driving (by their sheer number of particles), one would 
expect WR winds to be independent of their initial Z and the requirements of 
the collapsar model may never be met. If, on the other hand , iron (Fe) is largely 
responsible for the driving (as in O stars; Vink et al. 200 ll ). WR winds might 
indeed be less efficient in low Z environments. 



T o address the question regarding the Z dependence of WR winds. I Vink de Koter 
(120051 ) computed mass-loss rates for late-type WN and WC stars as a function 
of the initial metal content (representative of the host galaxy). The results are 
shown in Fig. [TJ For a discussion of the flattening in the mass-loss- Z dependence 
for initial metallicities below log (Z/Zq) = —2 an d potential c onsequences for 
the first stars (Pop III), the reader is referred to Vink ( 20061 ) . but for the Z 
range down to log (Z/Z & ) = —2, the mass loss is found to drop steeply, as M 
oc Z 0,85 , for the WN phase - where WR stars spend most of their time. This 
inefficiency of WR mass loss at subsolar Z may prevent the loss of stellar angular 
momentum, providing a boost to the collapsar model. 



3. Mass loss from OB stars: absolute rates and the bi-stability jump 



We switch from a discussion of mass-loss versus Z to a discussion of mass loss 
versus T c g. This is best described in terms of the wind efficiency number 
r\ = (Mwoo)/(L*/c), a measure for the momentum transfer from the photons 
to the ions in the wind. Here is the termin al velocity of the outflow and 



the luminosity of the star. IVink et al.l ((2000) computed wind models as a 



function of effective temperature as shown in Fig. [2j The overall behaviour is 
one of decreasing n with decreasing T e fj due to a growing mismatch between 
the wavelengths of the maximum opacity (in the UV) and the flux (moving to 
longer wavelengths). The behaviour reverses at the "bistability jump" (BSJ; e.g. 
Lame rs et al. 1995), wh ere r] increases by a factor of 2-3, as Feiv recombines to 
Fern (IVink et al.lll999T ). 

Recent mass- loss studies (|Trundle Lennonl 120051 : ICrowther etaT\ |2006) 
have reconfirmed d iscrepancies betw een empirical mass-loss rates and predictions 



for B su pergiants (Vink et al 



|2000| ). Discrepancies have also been reported for 



O stars (jFullerton et al.ll2006l ). and it is as yet unclear whether the reported 
discrepancies for B supergiants could be due to model assumptions (e.g. the 
neglect of wind clumping) or the physical reality of the BSJ. The most accurate 
way to derive M is believed to be through radio observations. Intriguingly, 
Benaglia et al. ( 20071 ) present empirical radio mass-loss rates as a function of 
effective temperature that show a similar behaviour to the mass-loss efficiency 
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Figure 2. Wind efficiency r\ = (Afz; oc ,)/(£ H ./ c) as a function of effective 
temperature. These predictions are taken from IVink et alj (|2000f ). Note the 
presence of the bistability jump around 25 kK, where r\ increases as Fe re- 
combines to Fein. 



predicted by IVink et al.1 (|2000l ). This may well be the first evidence of the 
presence of a mass-loss BSJ at the boundary between O and B supergiants. 
The relevance for stellar evolution is that when massive stars evolve to lower 
T e g after the O star main sequence phase, they are expected to cross the BSJ. 
Interestingly, LBVs brighter than log (L/L Q ) = 5.8 (see Fig. [3]) are expected to 
encounter it continuously - on timescales of their photometric variability, which 
we discuss in the next section. 



4. Mass loss from Luminous Blue Variables 



LBV s are unstable massive stars in the upper part of the HRD (e.g. lHumphrevs Sz Davidson 
1994). As can be seen in Fig. [3J the classical LBVs, like AG Car, are anticipated 
to cross the BSJ at ~ 21 000 K. One of the defining characteristics for LBVs is 
their S Doradus (SD) variation of ~1 — 2 mag on timescale s o f years (short SD 



phase s) to decades (long SD phases) ( van GenderenI 1200 ll ). IVink de Koter 
(2002) computed LBV mass- loss rates as a function of T P ff - shown in Fig. 
Overplotted are the empirical mass- loss rates for AG Car (|Stahl et a l. 2001), 
which vary on the times cales of the photometric variability. Although the agree- 
ment is not perfect (see I Vink Sz de Koterll2002l . for a discussion), the amplitude 
of the predicted variability fits the observations well, and most importantly the 
overall behaviour appears to be very similar, and may indeed be explained in 
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Figure 3. The LBVs in the HRD. The shaded areas represent the S Dor 
instability strip (diagonal) and the position of the LBVs during outburst 
(vertical). The dashe d vertical lin e at 21 000 K indicates the bistability jump. 
Figure adapted from lSmith et al.l (|2004f ) . 



terms of the physics of the BSJ . This bi-stable behaviour in an individual stellar 
wind (jPauldrach Pulslll990h causes the star to flip back and forth between 
two states: that of a low mass loss, high-velocity wind, to a high mass-loss, low 
velocity wind. The wind density (oc M/voo) would therefore be expected to 
change by a factor of ~2 x ~2, i.e. ~4 on the timescale of the SD variations. In 
the absence of any other material around the star, this would result in a pattern 
of concentric shells of varying density. 



5. Radio supernovae and progenitor mass loss 

Radio SNe (RSNe) lightcurves and the model for SN interaction w i th th e sur- 
rounding circumstellar material has been reviewed bv lWeiler et al. ( 1986), and 



is shown schematically in Fig. [5j The radio emission is due to non-thermal elec- 
trons, while the absorption may be due to both synchrotron self absorption as 
well as free-free absorption. Examples of the rise, peak, and power-law decline 
of radio lightcurves are shown in Fig. (The episodic bumps at late time are 
discussed in Sect. [6TJ) 

The model constrains the wind density and thu s the ratio of M to the ter- 
minal wind velocity: p oc M /v^r 2 . Assuming Voo, IWeiler et al. ( 20021 ) list M 



values in the range 10 -6 -10 -4 MQyr -1 . Fortunately, these values agree with 
mass-loss predictions, but are broadly representative for massive stars over al- 
most all post-main sequence evolutionary phases, making it hard to infer the 
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Figure 4. Predicted (dotted line) and empirical (dashed line) mass-loss 
rates versus T e g for the LBV AG Car. Note that both the qualitative be- 
havi our and the amp l itude of the mass-loss variations are well reproduced. 
See lVink fc de Koterl (pOOl for details. 



progenitor from radio lightcurves alone, unless these lightcurves betray their 
progenitor in some other way. 



6. Quasi-periodic oscillations in radio SNe lightcurves 

A number of recent RSNe have s hown sinusoidal m odulations in t heir radio 

lightc urves, in particular SN2001ig (|Rvder et al.ll2004l ) and SN2003bg (jSoderberg et al 
2006j) are strikingly similar in terms of both amplitude and variability timescale 
(see Fig. [6j). The recurrence tim escale t of the bumps is ~ 150 days. Using 
Eq. (13) from lWeiler et al.1 (|1986T ): 



AP 



-Rshcll 
^wind 



^ejecta ti ( t 
Wwind m \ti 



(1) 



where m is the deceleration parameter (here m = 0.85) and t{ is the time of 
measurement of the ejecta velocity relative to the moment of the explosion. 
Assuming v^md = 10-20 km sec -1 , typical wind velocities for red (super)giants, 
Ryder et al. ( 20041 ) found a period P between successive mass- loss phases that 
was too long for red (super)giant pulsations (100s of days), but too short for 
thermal pulses (10 2 -10 3 years). They therefore invoked an edge-on, eccentric 
binary scenario involving a WR-star and a massive companion. One of the main 
differences between LBV and red giant winds is that LBV winds are about 10 
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Figure 5. Model for SN ejecta interaction with the progenitor's wind. The 
radio emission arises at the interface between the outgoing shock and the most 
recent stellar wind. The various optical depth (r, r', t") contributions are 
respectively from a smo oth wind, a clumped wind and a potential intervening 
Hn region. Taken from iWeiler et all (|2002f) . 



times faster. If the progenitor of SN2001ig were an LBV, the expected period 
between successive mass-loss episodes would be AP ~ 2 5 yr (for an assume d 
^win d = 200 km sec -1 ) , cons istent with the long SD phase (jKotak Vinkll2006l ). 
Soderberg et al.l (|2006l ) infer density enhancements of a factor of ~2 dur- 



ing the deviations from pure power-law evolution. They consider a range of 
options that might account for the modulations, but they favour a single-star 
progenitor model of a WR star that underwent episodes of intensified mass loss. 
However, they do not specify the physical mechanism that gives rise to these 
periods of enhanced mass loss. Our SD mechanism for LBVs may alleviate this 
shortcoming. 



7. Discussion: do LBVs explode? 

Are LBVs viable SNe progenitors? It may be relevant that both SNe2001ig 
and 2003bg are "transitional" objects. SN2001ig was initially classified as type 
II (showing H lines) but metamorphosed into a type Ib/c object (no H lines, 
weak He lines) about 9 months later. This suggests that it has lost most of its 
H-rich envelope. SN2003bg however was first classified as a type Ic, but within 
a month the spectrum evolved into a type II SN. This transitional behaviour 
hints at the fact that their progenitors are intermediate evolutionary objects: 
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Figure 6. Radio luminosity versus time for two strikingly similar recent 
SNe: 2001ig and 2003bg. Note the quasi-sinusoidal modula tions during the 
power-law decline phase. Taken from lSoderberg et al.l (|2006[ ). 



H-rich compared to OB/red (super)giants, but H-poor compared to WR stars. 
LBVs are likely candidates. 

During this meeting, there was discussion about clumping in O star winds. 
The value for the clumping factor in O star winds is very much an open issue. If 
these factors would be much larger than two, this would have severe implications 
for massive star evolution. One consequence might be that giant LBV eruptions 
(77 Car type eruptions, not t he typifying SD variations) dominate the integrated 
mass loss during evolution ( Smith! 20061 ). An alternative scenario could be that 
post-main sequence stars do not become WR stars, but explode early - during 
their LBV phase. 

Here, we have presented indications that at least those SNe that show quasi- 
periodic modulation s in their radio lightcurves might have LBV progenitors 
( Kotak Sz Vink 2006). It has also been speculated that LBVs may be the generic 
progenitors of type Iln SNe (?), however the type Iln phenomenon (arising from 
SN ejecta expanding into dense circu mstellar matter) m ay be relevant to both 
core-collapse and thermonuclear SNe ( Kotak et al.ll2003 ). and although it may 
be reasonable to expect that some type Iln SNe have LBV progenitors, there 
remains a lot of work to be done to prove this. Nonetheless, it appears that 
the "standard scenario" for massive star evolution may need revision. Future 
mass-loss predictions will certainly play a major role in adjusting even our most 
basic knowledge of massive star evolution. 
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